On the coherent rotation of diffuse matter in numerical simulations of galaxy clusters

We present a study on the coherent rotation of the intracluster medium and dark matter components of simulated galaxy clusters extracted from a volume-limited sample of the MUSIC project. The set is re-simulated with three different recipes for the gas physics: $(i)$ non-radiative, $(ii)$ radiative without AGN feedback, and $(iii)$ radiative with AGN feedback. Our analysis is based on the 146 most massive clusters identified as relaxed, 57 per cent of the total sample. We classify these objects as rotating and non-rotating according to the gas spin parameter, a quantity that can be related to cluster observations. We find that 4 per cent of the relaxed sample is rotating according to our criterion. By looking at the radial profiles of their specific angular momentum vector, we find that the solid body model is not a suitable description of rotational motions. The radial profiles of the velocity of the dark matter show a prevalence of the random velocity dispersion. Instead, the intracluster medium profiles are characterized by a comparable contribution from the tangential velocity and the dispersion. In general, the dark matter component dominates the dynamics of the clusters, as suggested by the correlation between its angular momentum and the gas one, and by the lack of relevant differences among the three sets of simulations.Comment: 12 pages, updated to match the MNRAS versio

Morphological estimators on Sunyaev-Zel'dovich maps of MUSIC clusters of galaxies

The determination of the morphology of galaxy clusters has important repercussions for cosmological and astrophysical studies of them. In this paper, we address the morphological characterization of synthetic maps of the Sunyaev-Zel'dovich (SZ) effect for a sample of 258 massive clusters (Mvir> 5×1014h-1M⊙at z=0), extracted from theMUSIC hydrodynamical simulations. Specifically, we use five known morphological parameters (which are already used in X-ray) and two newly introduced ones, and we combine them in a single parameter. We analyse two sets of simulations obtained with different prescriptions of the gas physics (non-radiative and with cooling, star formation and stellar feedback) at four red shifts between 0.43 and 0.82. For each parameter, we test its stability and efficiency in discriminating the true cluster dynamical state, measured by theoretical indicators. The combined parameter is more efficient at discriminating between relaxed and disturbed clusters. This parameter had a mild correlation with the hydrostatic mass (~0.3) and a strong correlation (~0.8) with the offset between the SZ centroid and the cluster centre of mass. The latter quantity is, thus, the most accessible and efficient indicator of the dynamical state for SZ studiesThis work has been partially supported by funding from Sapienza University of Rome - Progetti di Ricerca Anno 2015 prot. C26A15LXNR. GY and FS acknowledge financial support from MINECO/FEDER under research grant AYA2015-63810-P. ER acknowledge financial contribution from the agreement ASI-INAF n 2017-14-H.

Machine learning methods to estimate observational properties of galaxy clusters in large volume cosmological N-body simulations

This is a pre-copyedited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. The version of record Monthly Notices of the Royal Astronomical Society 518.1 (2023): 111-129 is available online at: https://academic.oup.com/mnras/article-abstract/518/1/111/6795309?redirectedFrom=fulltext#no-access-messageIn this paper, we study the applicability of a set of supervised machine learning (ML) models specifically trained to infer observed related properties of the baryonic component (stars and gas) from a set of features of dark matter (DM)-only cluster-size haloes. The training set is built from the three hundred project that consists of a series of zoomed hydrodynamical simulations of cluster-size regions extracted from the 1 Gpc volume MultiDark DM-only simulation (MDPL2). We use as target variables a set of baryonic properties for the intracluster gas and stars derived from the hydrodynamical simulations and correlate them with the properties of the DM haloes from the MDPL2 N-body simulation. The different ML models are trained from this data base and subsequently used to infer the same baryonic properties for the whole range of cluster-size haloes identified in the MDPL2. We also test the robustness of the predictions of the models against mass resolution of the DM haloes and conclude that their inferred baryonic properties are rather insensitive to their DM properties that are resolved with almost an order of magnitude smaller number of particles. We conclude that the ML models presented in this paper can be used as an accurate and computationally efficient tool for populating cluster-size haloes with observational related baryonic properties in large volume N-body simulations making them more valuable for comparison with full sky galaxy cluster surveys at different wavelengths. We make the best ML trained model publicly availableThe authors thank the anonymous referee for his/her invaluable comments and suggestions, without which this work would be in complete. D.d.A., W.C. and G.Y. would like to thank Ministerio de Ciencia e Innovación for financial support under project grant PID2021-122603NB-C21. WC is supported by the STFC AGP Grant ST/V000594/1 and the Atracción de Talento Contract no. 2020-T1/TIC-19882 granted by the Comunidad de Madrid in Spain. He also thanks the Ministerio de Ciencia e Innovación (Spain) for financial support under Project grant PID2021-122603NB C21. He further acknowledges the science research grants from the China Manned Space Project with NO. CMS-CSST-2021-A01 and CMS-CSST-2021-B01. G.M. acknowledges financial support from PID2019-106827GB-I00/AEI / 10.13039/501100011033 The CosmoSim database used in this paper is a service by the Leibniz Institute for Astrophysics Potsdam (AIP). The MultiDark database was developed in cooperation with the Spanish MultiDark Con solider Project CSD2009-00064. The authors acknowledge The Red Española de Supercomputación for granting computing time for running the hydrodynamical simulations of The300 galaxy cluster project in the Marenostrum supercomputer at the Barcelona Super-computing Cente

nIFTy galaxy cluster simulations – II. Radiative models

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $\Lambda$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction.Comment: 20 pages, 13 figures, submitted to MNRA

The MUSIC of Galaxy Clusters I: Baryon properties and Scaling Relations of the thermal Sunyaev-Zel'dovich Effect

We introduce the Marenostrum-MultiDark SImulations of galaxy Clusters (MUSIC) Dataset, one of the largest sample of hydrodynamically simulated galaxy clusters with more than 500 clusters and 2000 groups. The objects have been selected from two large N-body simulations and have been resimulated at high resolution using SPH together with relevant physical processes (cooling, UV photoionization, star formation and different feedback processes). We focus on the analysis of the baryon content (gas and star) of clusters in the MUSIC dataset both as a function of aperture radius and redshift. The results from our simulations are compared with the most recent observational estimates of the gas fraction in galaxy clusters at different overdensity radii. When the effects of cooling and stellar feedbacks are included, the MUSIC clusters show a good agreement with the most recent observed gas fractions quoted in the literature. A clear dependence of the gas fractions with the total cluster mass is also evident. The impact of the aperture radius choice, when comparing integrated quantities at different redshifts, is tested: the standard definition of radius at a fixed overdensity with respect to critical density is compared with a definition based on the redshift dependent overdensity with respect to background density. We also present a detailed analysis of the scaling relations of the thermal SZ (Sunyaev Zel'dovich) Effect derived from MUSIC clusters. The integrated SZ brightness, Y, is related to the cluster total mass, M, as well as, the M-Y counterpart, more suitable for observational applications. Both laws are consistent with predictions from the self-similar model, showing a very low scatter. The effects of the gas fraction on the Y-M scaling and the presence of a possible redshift dependence on the Y-M scaling relation are also explored.Comment: 22 pages, 25 figures, accepted for pubblication by MNRA

nIFTy galaxy cluster simulations - IV. Quantifying the influence of baryons on halo properties

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern smoothed particle hydrodynamics (SPH) approaches. For each code represented we have a dark-matteronly (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<≈ 20 per cent) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include active galactic nucleus feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers

nIFTy galaxy cluster simulations – I. Dark matter and non-radiative models

We have simulated the formation of a galaxy cluster in a Ʌ cold dark matter universe using 13 different codes modelling only gravity and non-radiative hydrodynamics (RAMSES, ART, AREPO, HYDRA and nine incarnations of GADGET). This range of codes includes particle-based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span classic and modern smoothed particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at z = 0, global properties such as mass and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes RAMSES, ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing classic SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid-based methods

Large scale numerical simulations of the Sunyaev-Zel’dovich effect: the next generation tool for Precision Cosmology

Tesis doctoral inédita leída en la Universidad Autónoma de Madrid, Facultad de Ciencias, Departamento de Física Teórica. Fecha de lectura: 27-05-201

Baryon and Sunyaev-Zel’dovich Effect properties of MareNostrum and MultiDark Simulated Clusters (MUSIC)

We report the first results of the MUSIC project. It consists of two data sets of resimulates clusters extracted from two large dark matter only simulations: Marenostrum Universe and Multidark. In total, the MUSIC contains more than 400 clusters resimulated with high resolution both with radiative and non-radiative physics included. Here we present the first results on the properties of the baryon content and the Sunyaev Zeldovich scaling relations

The evolution of the Y-M scaling relation in MUSIC clusters

This work describes the baryon content and Sunyaev-Zeldovich properties of the MUSIC (Marenostrum-MultiDark SImulations of galaxy clusters) dataset and their evolution with redshift and aperture radius. The MUSIC dataset is one of the largest samples of hydrodynamically simulated galaxy clusters (more than 2000 objects, including more than 500 clusters). We show that when the effects of cooling and stellar feedbacks are properly taken into account, the gas fraction of the MUSIC clusters consistently agrees with recent observational results. Moreover, the gas fraction has a net dependence with the total mass of the cluster and increases slightly with redshift at high overdensities. The study of the Y-M relation confirms the consistence of the self-similar model, showing no evolution with redshift at low overdensities. (C) 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinhei